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Quantum dynamics of a two-state system induced by a chirped zero-area pulse 
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It is well known that area pulses make Rabi oscillation and chirped pulses in the adiabatic in¬ 
teraction regime induce complete population inversion of a two-state system. Here we show that 
chirped zero-area pulses could engineer an interplay between the adiabatic evolution and Rabi-like 
oscillations. In a proof-of-principle experiment utilizing spectral chirping of femtosecond laser pulses 
with a resonant spectral hole, we demonstrate that the chirped zero-area pulses could induce, for 
example, complete population inversion and return of the cold rubidium atom two-state system. 
Experimental result agrees well with the theoretically considered overall dynamics, which could be 
approximately modeled to a Ramsey-like three-pulse interaction, where the x- and z- rotations are 
respectively driven by the hole and the main pulse. 

PACS numbers: 32.80.Qk, 32.80.Wr, 42.65.Re 


I. INTRODUCTION 

Quantum technology, a new field of engineering based 
on quantum mechanics, is expected to make a defining 
impact on life in the 21st century [lj. In many appli¬ 
cations of quantum technology, such as quantum com¬ 
puting, quantum cryptography, quantum simulation, and 
quantum metrology, the ability to control the dynamics 
of quantum systems is an essential and foremost impor¬ 
tant necessity. Often referred to as a qubit in quantum 
information science [2], which is used as an elementary 
building block in more complex quantum devices and ma¬ 
chines, a two-state quantum system is defined by a co¬ 
herent superposition of two energy states. 

The dynamics of a two-state quantum system induced 
by coherent radiation have been extensively studied: 
Rabi oscillations m, Ramsey fringes [5j , and spin-echo 
dynamics [Bj are among the best known examples. These 
examples are well understood in the basis of the area 
pulse concept 0 0, where the pulse (temporal) area 
0 is defined by 0 = f^° f l(t)dt, with Rabi frequency 
17 (t) that is defined (e.g., in laser-atom interactions) by 
17(t) = p,£(t)/H with /i the atomic transition dipole mo¬ 
ment and £ (t) the envelope of a laser electric-field in res¬ 
onant frequency. Besides Rabi oscillations, there are only 
a few exactly solvable models, including Landau-Zener, 
Rosen-Zener, Allen-Eberly, Bambini-Berman, Demkov, 
Nikitin, and Carrol-Hioe models (4j I9HT2] . 

In general, coherent radiation is a complex, therefore 
powerful, control means; its amplitude, phase, frequency, 
and polarization can be used as independent control pa¬ 
rameters. Of particular relevance in the context of the 
present paper, another important control parameter in 
two-state system dynamics is chirp mm, the frequency 
sweeping rate of coherent radiation in time. When a laser 
frequency co(t) changes linearly in time with the (tempo¬ 
ral) chirp parameter a as u>(t) = w 0 + 2 at, across w D , 


’Electronic address: jwahnOkaist.ac.kr 


the Bohr transition frequency of the two-state system, 
a complete population inversion (CPI) between the two 
states occurs (when the adiabatic condition [15] is satis¬ 
fied all along the interaction), and the two-state system 
evolves in time through a so-called rapid adiabatic pas¬ 
sage (RAP) [I5j . 

In this paper, we consider chirped zero-area pulses to 
study the two-state system dynamics induced by them. 
To begin with, in the basis of the pulse-area theorem m, 
“simple” zero-area pulses make no net transition, leaving 
the two-state system intact. However, to be more precise, 
the system does change during the dynamics, although 
it ends up as a complete population return (CPR) after 
the completion of the interaction. In that regards, cer¬ 
tain manipulation of the initial zero-area pulse could alter 
the system evolution completely, leading to a significant 
change in the final transition probability. Even CPI could 
take place by detuning the zero-area pulse mm, where 
abrupt phase change makes transitions in adiabatic basis 
or a crossover between adiabatic evolution and Rabi os¬ 
cillations. Inspired by these counter-intuitive examples, 
we proceed further to consider chirping an (on-resonant) 
zero-area pulse in this paper. For this, we have con¬ 
ducted laser-atom interaction experiments with intense 
shaped laser pulses and spatially confined atomic vapor. 
The results are (1) chirped zero-area pulses could produce 
both CPI and CPR; (2) the excited-state probability os¬ 
cillates as a function of the effective pulse-area defined 
by pulses with zero chirp; and (3) the given dynamics 
can be modeled in terms of Ramsey-type three-pulse se¬ 
quence R x (tt/2)R z (Q z )R x (tt/2) (to be explained below), 
all of which can be summarized by an interplay between 
adiabatic evolution and Rabi-like oscillations. 

The rest of the paper is organized as follows: In Sec. 
II, we describe the model and theoretically analyze the 
coherent excitation induced by chirped zero-area pulses, 
and the detailed interpretation of the dynamics will be 
considered in Sec. III. Section IV is devoted to the ex¬ 
perimental description, and we present the experimental 
results in Sec. V, before concluding in Sec. VI. 
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II. THEORETICAL CONSIDERATION 

When the electric field of a laser pulse is defined in 
time domain as 

E(t) = £{t) cos (u} 0 t + <p), (1) 

where £(t) and p are the envelope and phase, respec¬ 
tively, of the electric-field that oscillates with the carrier 
frequency of w D set to the Bohr transition frequency of 
the two-state atom, the pulse area 0 is given under the 
rotating wave approximation by 

B = \ r E(t)e~ iUot dt = (2) 

n J-oo n 

where E{u>) is the amplitude spectrum of the electric 
field. In that regards, on-resonance pulses with zero 
resonant-frequency component, i.e., E{oj 0 ) = 0, are zero- 
area pulses (0 = 0 ) defined in frequency domain, which is 
consistent with the zero-area pulse alternatively defined 
in time domain m- For instance, a Gaussian pulse with 
a spectral hole around the resonance frequency oj 0 is a 
zero-area pulse, which is given by the difference of two 
Gaussian pulses both frequency-centered at ui a , i.e., 


with r o i = 2/Auj.i, the initial pulse widths, and t* = 
T oi + 4 c |/ t 2 the chirped pulse widths. The chirp pa¬ 
rameters are a = 2c2/ / (tq 1 -|-4c|) and f) = 2c 2 /(tq j2 +4c2), 
respectively, for the main and hole pulses. Note that, al¬ 
though Eq. <[5j) satisfies 0 = 0 based on the definition 
in Eq. ([2]), the conventional pulse area defined in time 
domain,/x£ to ta \{t)/Mt, for real envelope £total (t) of 
the total electric field E(t) = Ei{t) — E 2 {t), is generally 
non-zero for c 2 7 ^ 0 . 

As to be discussed below, this type of chirped zero- 
area pulses can induce both CPI and CPR of a two-state 
system. To understand the nature of the given dynam¬ 
ics, we will transform the two-state system in three steps: 
first to an interaction basis (to freeze the phase evolution 
of the bare atomic state), then to an adiabatic basis [T5] 
(to trace the system’s state-mixing during the interac¬ 
tion), and finally to another interaction basis (to clear 
out the remaining fast phase-evolution of the state). Af¬ 
ter obtaining the model Hamiltonian, we will discuss the 
resulting state-evolution in terms of Ramsey-like three 
rotations in the Bloch sphere representation. 

A. Schrodinger equation described in the 
interaction picture of the main pulse 


E (v) 

E n 


= exp 


(w - w 0 ) 2 ' 
Aw 2 . 


— exp 


(a) - UJq) 2 ' 

A - 


(3) 


When the spectral width of the first pulse, Awi, is sig¬ 
nificantly bigger than that of the second, Aoj 2 , {i.e., 
Awi A lo 2 ), the second pulse (or, to say, the hole pulse) 
is regarded as a narrow spectral hole of the first pulse (the 
main pulse). 

Let us consider now a chirped zero-area pulse, by chirp¬ 
ing the above zero-area pulse, which reads 


E(u) 


■ 

(w-w 0 ) 2 


(w-w 0 ) 2 T 

exp 

Aw 2 

— exp 

A 10 % \ _ 


x exp 



(4) 


The dynamics of the two-state system that interacts 
with the electric field E(t) in Eq. ([ 5 ]) is governed by the 
interaction Hamiltonian Vj as 

0 n 2 (t)' 

n* 2 {t) 0 

= V Itl + V L2 , (8) 

where A (t) = uj 0 — ui(t) = —2 at is the (instantaneous) 
detuning EDI- The Rabi frequencies of the main and hole 
pulses are, respectively, given by 


£li{t) — 

fat), 

(9) 

fl 2 {t) = 

-£ 2 (i)e i[(/3 ~ a)t2_M , 

h 

(10) 


h 

—A(t) 

h 

2 

fli(t) A(t) 

~ 2 


where the two constituent Gaussian pulses are simulta¬ 
neously chirped with the frequency chirp c 2 . The electric 
field in time domain is then given by 


where we define = <p\ — p> 2 , and the time-independent 
phase tp\ is included in the base vector for the sake of 
simplicity. 


E(t) = £i(t) cos + at)t + p>\] 

-£ 2 (0 cos [{u 0 + /3t)t + ip 2 ] 

= E^t) - E 2 {t), (5) 

where the amplitude and the phase for each pulse * = 1,2 
are, respectively, 


£i{t) = 


/ T 0) j t 2 / 2 

(6) 

V2 ) 

n 

Vi = 

1 

— tan 

2 

-i 2 c 2 

0,1 

(7) 


B. Dynamics described in the adiabatic basis of 
the main pulse 

In Eq. ([8]), Vqi, the first part of the interaction Hamil¬ 
tonian, varies slowly compared to the second Vi : 2 , so the 
system dynamics can be more easily understood in the 
adiabatic basis [1]J] of the main pulse. The eigenvalues 
of Vip are given by 

h ± {t) = ±lsJni{t) + A\t) 


(ii) 
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and the corresponding eigenstates are 

| = cosi?(£)|0)j — sini?(£)|l)j, 

| ip+{t)) = sini?(t)|0)j+cosi?(t)|l)/, (12) 

where |0)/ and |l)j form the eigenbasis in the interaction 
picture, and the mixing angle d[t) is defined by 


the phase factor, exp(i f* \A(t')\dt') = exp(z|A(£)|£/2), 
through a transformation to another interaction basis, 
i.e., 

mt)) F = TMt)) A , ( 20 ) 

where Ta = exp (i f* H A (t')dt'/hj with 


d{t) = ^ tan 1 for 0 < d(t) < (13) 


H a = 


h 

2 


— |A(t)| 0 

0 | A(t) | 


( 21 ) 


So, the two-state system can be described in adiabatic ba¬ 
sis (IV’-(t)), \4>+(t))) with the transformation |^>(£)}a = 
R{fl)\ip(t))i, where the rotation i?($) is defined by 


m 


cos $(£) — sind(t) 
sind(t) cos d(t) 


(14) 


Note that each adiabatic base vector, \i/j_(t)) or |^> + (£)), 
changes from one atomic state to the other, as time 
evolves from t = —oo to oo. In other words, when C 2 > 0, 


lim |= |0)/, lim | ip-(t)} = -|1)/, 

t —>—oo £—> oo 

lim | ip+(t)) = |l)j, lim | ip+{t)) = |0)/, (15) 

£—>■ — 00 £—> oo 


and when C 2 < 0, the relation is reversed, since the time 
dependence of A is opposite. 

In the given adiabatic basis, the Schrodinger equation 
is given by 

ih^\ip(t))A = (RViFr 1 + iHRR~ 1 ^ | i>(t))A, (16) 


where the second term in the parenthesis can be ignored, 
when the evolution by S2i(f) is adiabatic, i.e., 


ihRR 1 = ih 


0 -i9 

d 0 


0, 


(17) 


so the resulting interaction Hamiltonian is given by 


V A 


h [ a_ -n 2 

2 -fi 2 A + 


+ Mi(fl 2 ) sind 


where 3fJ(fl 2 ) is the real part of f l 2 (t). 


cosd 

sind 


sind 

— cos d ’ 


(18) 


C. Ramsey-type three pulsed interactions 


Let us take a closer look at Hamiltonian Va- At the 
extreme times |£| —> oo, it becomes 


V A = a 


—IA(t)| -n 2 (t) 
-n* 2 (t) |A(t)| 

-|A(t)| nm 
n 2 (t) |A(t)| 


for t —> —oo 


for t —>• oo. (19) 


Note that the phase evolutions of the diagonal terms are 
opposite with each other. So, it is convenient to remove 


The resulting final Hamiltonian becomes 


V F = ?a(Va - H A )Tl = ^ 


—A F {t) Sl F (i) 

A F (t) \ ’ 


( 22 ) 


where the effective “detuning” A F and “coupling” Op- 
are given by 

A f = \J O \{t) + A 2 (t) + 5J(0 2 ) sin[2d(£)] - |A(£)|, 
0 F = {-cos[2d(t)]3fi(02)-i3(02)}e- l|A(t)|t/2 , (23) 


where 3(02) is the real part of fl 2 {t). The above Hamil¬ 
tonian leads to the detuned Rabi oscillation m after 
a second rotating-wave approximation. The remaining 
time-dependence in 0 F (f), the slowly-varying phase fac¬ 
tor is negligible for the considered chirp values. 

Figure ITJa) plots the Rabi frequencies of the main 
pulse, Oi(f), the hole, |0 2 (£)|, and the total electric field, 
|0(t)| = |Or(£) — 0 2 (£)|, of a chirped zero-area pulse that 
results in a CPR. The numerical values for the given 
pulse are c 2 = 5.1 x 10 4 fs 2 , Aoq = 1.5 x 10 13 rad/s, 
and Aa> 2 = 1.9 x 10 12 rad/s. The peak electric field of 
the main pulse is E 0 Au)\/y/2 = 2.3 x 10 8 V/m, and the 
dipole moment of atomic rubidium ( 85 Rb) for linearly 
polarized light is given by /i = 1.46 x 10~ 29 Cm Tl 1 3i?l . 
To examine adiabaticity, we use the adiabaticity function 
/(£) = |S2(£)A(£) — fl(£)A(£)|/2[\/S2 2 + A 2 ] 3 [15] . Adia¬ 
batic evolution condition is given by f(t) -C 1. For lin¬ 
early chirped pulses, the function /(f) can be explicitly 
written as, 


h(t) 

h{t) 


q|fli(f)|(2f 2 /T 1 2 + 1) 

[v/1STp| 2 + 4a 2 £ 2 ] 3 ’ 
/3|H 2 (f)|(2t 2 /r 2 + l) 
[v/|H 2 | 2 +4/3 2 f 2 ] 3 ’ 


(24) 

(25) 


where /i(t) and f 2 (t) are the adiabatic functions for the 
main and the hole pulses, respectively. Note that the 
pulse width r,; and the chirp parameters a and /3 are 
functions of the bandwidth Axy as well as of the linear 
chirp C 2 - With the given parameters in Fig. [TJ we obtain 
the adiabatic condition, /i(£) < 0.2, for the main pulse 
with E 0 Au>i/y/2 > 1.2 x 10 8 V/m and C 2 > 2 x 10 4 fs 2 , 
which ensures the adiabatic evolution by the main pulse. 
However, we obtain f 2 {t) > 1.2, so the hole pulse induces 
simple Rabi oscillations, even up to the region C 2 — 1 x 
10 5 fs 2 . 
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In Fig. 0 b), the time-evolution of the transition prob¬ 
ability in adiabatic basis is compared with the transition 
in the bare basis, which clearly shows the transition in 
the adiabatic basis, or a CPR. Also, we plot in Fig. [ljc) 
the polar angle, 0 lot (t ), of the rotational axis on the Bloch 
sphere in the adiabatic basis, defined by 

0 Iot (t) = tan" 1 -L^jj. (26) 

The rotation on a Bloch sphere is approximately a phase 
evolution (about the z axis) when 9 rot = 0, and Rabi 
oscillations when 0 Tot = 0.57T. Therefore, the plateaus in 
Fig. [ijb) suggest that there exist three distinct coupling 
regimes: In the two tail regions ( t < —t\ and t > ti), 
the “coupling” is dominant, so a Rabi oscillation is ex¬ 
pected; and, in the central region (— t\ < t < Ti), the 
“detuning” is large, effectively causing approximately a 
phase evolution of the system. 


1. |A(t)| Sli(t) (the tail regions) 


In the tail regions (t < — n and t > n), the detuning 
A(f) greatly exceeds the main-pulse interaction fli(i), 
i.e., |A(f)| fli(t). So, we get Ap —> 0, and f If —► 
-|fl 2 |e^ (for t<n) or (for t > n). Given 

that, the Hamiltonian in Eq. \22\ is approximated, when 
the slowly varying phase e i/3t in Q 2 is neglected, as 


Vf 


0 | n 2 {t)\e i( - n+ ^ 

\Vl 2 {t)\e- i ^ + ^ 0 


for t < — n 


0 

\n 2 (t)\t 


|H 2 (t)|t 

o 


for t > T \, 


(27) 


which means that the dynamics in the tail regions are 
rotations on the Bloch-sphere surface, represented by 
R n p v (@ 2) and R_ v (Qt), respectively, for t < t\ and 
t > t-[. The rotation axes are on the xy plane, of which 
the directions are defined by azimuthal angles <j> = tt + if 
and 4> = —if, respectively. The rotation angles are given 

by 


/ —Tl POO 

\n 2 (t)\dt and 0^" = / \Q 2 (t)\dt, (28) 

-OO J Tl 

as a function of the hole pulse E 2 (t), and 0^=0^" due 
to the symmetry. 


2. |A(t)| fli (t) (the central region) 

In the central time region (—Ti < t < Ti), the main 
pulse Ei(t) plays an important role. To understand 
the dynamics in this region, we consider an extreme 
approximation under the conditions A(f) and 

H^t) |fl 2 (t)|. (Note that the case depicted in Fig. 
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FIG. 1: (Color online) (a) The Rabi frequencies fb(f) (dotted 
line) for the main pulse, | ST2 (t) | (dashed line) for the hole, 
and |fl(f)| (solid line) for the total electric field, (b) The 
time evolution of the transition probabilities in bare atom 
and adiabatic bases, (c) The transience of the polar angle 
8 rot (t) (solid line) for the rotational axis of the Bloch vector 
in the adiabatic basis compared with |A.F(t)| (dotted line) 
and | f (t)| (dashed line). The calculation parameters chosen 
at a CPR are given in the text. 


1 is different from this extreme case.) Then, we get 
A f ~ Hi(f) — |A(t)|, which leads to 


h r —Or(t) + |A(t)| 0 

2 0 fii(i)-|A(t)| 


(29) 


The resulting two-state dynamics is the rotation about 
the 2-axis, or R z ((-)-[), by an angle 

0!= r -\A(t)\)dt (30) 

J — Tl 

defined by the main pulse E\(t). 


As a result, the overall dynamics of the two-state sys¬ 
tem by the chirped zero-area pulse (in the transformed 
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adiabatic basis) can be summarized as following: 

i?_^(©+)i? z (0i)i? w+v (0^)., (31) 

which is not in the bare basis. This manifestation of the 
three-step interaction is a reminiscence of the Ramsey ro¬ 
tation often described by R x (tt/2)R z (Q z )R x (tt/2), from 
which the difference is to be discussed. 

Note that the necessary condition required (in our con¬ 
sideration) for Ramsey-type three pulse sequence (First 
Rabi + Second adiabatic evolution + Third Rabi) mod¬ 
eling is in fact the contrast in the spectral width between 
the main and hole pulses. When the spectral width of the 
main pulse is large enough, chirping makes the main pulse 
satisfy the adiabatic condition, while the same amount of 
chirp makes the spectrally-narrow hole pulse remain as 
a Rabi-inducing pulse. Therefore, pulses of smooth en¬ 
velope shapes other than Gaussian, when they have the 
spectral hole and the chirp, could also be approximated 
to a Ramsey-type three pulse sequence. 


III. CALCULATION RESULTS 


Figure [2] presents the numerical simulation for the two- 
state system dynamics induced by the chirped zero-area 
pulses. For the numerical calculation, the Schrodinger 


equation in Eq. (16) is used with Ao>i = 1.5 x 10 13 rad/s 


and Aw 2 = 1.9 x 10 i2 rad/s. In Fig. [2])a), the probability 
of the transition to the excited state in bare atomic ba¬ 
sis is shown as a function of 0R a bi and C 2 , where 0R a bi 
denotes the effective Rabi oscillation phase for the main 
pulse with zero chirp, defined by 


0R a bi = t I £i(t;c 2 = 0 )dt. 

n J-00 


(32) 



0 n 2n 3n 4n 5n 

©Rabi 




Note that 0R a bi was varied by changing E a , the electric 
field amplitude, in experiment and calculation. 

We can identify several CPR and CPI regions in Fig. [2] 
The CPR regions appear as localized dark spots, be¬ 
sides the wide dark region in low pulse energies; and the 
CPI regions are bright strips. To understand the dy¬ 
namics, we select a few characteristic points in Fig. [2](a) 
and trace their quantum trajectories in adiabatic basis. 
In Figs. §b, c), the trajectories are plotted on the ( 0 , 0 ) 
surface of the Bloch sphere (defined in adiabatic basis), 
where 0 and 0 denote the polar and azimuthal angles, 
respectively, of the Bloch vector. At the CPR points, 
marked by @, ©, and ©, each trajectory starts from 
0 = 0 and ends at 0 = n, as shown in Fig. [2](b). Like¬ 
wise, at the CPI points, marked by ®, (2), • • •, and ©, 
each trajectory starts from 0 = 0 and ends at 0 = 0 , as 
shown in Fig. [2^c). Note that CPR in bare atomic ba¬ 
sis, |0) -A |0), is CPI in adiabatic basis, |0_(— 00 )) —>■ 
|0 + (oo)), and CPI in bare atomic basis, |0) -A — 11), is 
CPR in adiabatic basis, |0_(— 00 )) -A |0_(oo)). 

In particular, each CPR trajectory (CPI in adiabatic 
basis) in Fig. [2jb) consists of three distinctive rotations: 


FIG. 2: (Color online) (a) Numerical calculation of the tran¬ 
sition probability in the atomic basis is plotted as a function 
of ©Rabi, the Rabi phase of pulses with zero chirp, and C 2 , 
the frequency chirp. Localized CPR regions appear as spots 
around @, (B), • • ■, (E); and CPI regions as strips, e.g., along 
(T)-(B). The dotted lines indicate contours for the estimated 
rotation angles Oi and ©2 defined by Eqs. (28) and (301, 
(b) and (c) Trajectories of CPR and CPI represented in the 
adiabatic basis for chosen points in (a), respectively. 


a rotation of A0 = 7t/2 about a rotation axis in the xy- 
plane, a wobbling rotation about the z-axis, and again a 
A0 = 7 r /2 rotation about another axis in the xy-plane. 
Note that the azimuthal rotation angles for the second 
wobbling rotations for the @, ®, © trajectories are 
about A0 = 27 t, 47 t, and 67 r, respectively. The CPR 
trajectories are therefore represented in accordance with 
Eq. ([31]) by 


R- lp (ir/2)R z (2nir + 'y)R n+(p (ir/2), 


(33) 


where n is a positive integer. It is noted that 7 = 
—if — (n + tp) = 7r — 2ip compensates the angle differ- 
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ence between the first and third rotation axes, explaining 
the deviation of the azimuthal rotation angle from 2mr. 
Here, from Eqs. ([T]) and (10), y>is given by ip = ipi — ip? = 
[tan^ 1 ( 2 c 2 /rQ a ) —tan^ 1 ( 2 c 2 /rg 2 )]/ 2 . Since Awi > Au; 2 , 
7 varies from 7r to ~7r/2, while c 2 varies from 0 fs 2 to 
100000 fs 2 , and 7 ~ 0 . 67 T at @, ®, and <Q. Similarly, 
the trajectories for © and © (not shown) can be approxi¬ 
mately represented by f?_ ¥ ,(37r/2)f? z (27r + 7 )f? w +¥>(37 t/ 2) 
and R_ v (3n/2)R z (4n+'y)R- K + lp {3n/2), respectively. 

The actual trajectories in Fig. §*>) are slightly de¬ 
viated from the approximated model in Eq. ©. It 
can be explained that, in the central region, the rota¬ 
tion axis is not perfectly aligned with the 2 -axis, caus¬ 
ing the wobbling of the 2 -rotation. As a result, CPRs 
do not occur at exact ©2 = mr + tt/2. Also, in the 
tail regions, the mixing angle d(f) does not converge to 
0 (for t < —n) or 7 r /2 (for t > 17 ) during the evolu¬ 
tion, and gives cos[2d(f)] = 1 — 5 (for t < —n) and 
cos[2$(i)] = — 1 + S (for t > n) in Eq. ( [23] ) , where S > 0 
is a small deviation. Since the azimuthal angle of a ro¬ 
tation axis is given by arg[f 2 ^(t)], the actual azimuthal 
angle difference between the two rotation axes becomes 
smaller than the model. Thus, as shown in Fig. H*), 
the CPRs (@, @, and (Q) occur near ©i = 2nn + 0 . 57 T, 
rather than ©! = 2mr + 7 (~ 2nn + 0 . 67 t). 

The intermediate region between the tail regions and 
the central region, appearing in Figs. [2]b,c) around the 
end the vertical evolution and beginning of the wobbling 
2 -rotation, is short in time compared to other regions (the 
ratio is roughly 0.2 compared to the central region). Also, 
the change of the rotation axis is much faster (roughly 
5 times faster) than the rotation. Therefore, as shown 
in Figs. |2]b,c), there are no significant influence to the 
dynamics in the intermediate region. 


IV. EXPERIMENTAL SETUP 

The experimental investigation of the chirped-zero 
area-pulse interaction of the two-state system is per¬ 
formed with shaped intense laser pulses and cold atomic 
vapor. The experimental setup is shown in Fig. |3] Femto¬ 
second laser pulses were first produced from a home¬ 
made Ti:sapphire laser oscillator with 80 MHz repeti¬ 
tion rate, and amplified about 5 x 10 5 times by a home¬ 
made Ti:sapphire multi-pass amplifier operating at 1 kHz 
repetition rate. The wavelength of the pulses was cen¬ 
tered at X 0 = 794.7 nrn, resonant to the atomic rubid¬ 
ium ( 85 Rb) 5 Si / 2 —> 5 Pi /2 transition. The two states 
of the quantum system were | 0 ) = |5Si/ 2 ) and | 1 ) = 
|5P 1/2 ). The laser bandwidth was AAfwhm = 7 nrn 
(Acc-| = 1.5 x 10 13 rad/s), which was narrow enough to 
restrict the quantum system to a two-state system, and 
the corresponding pulse duration was r 01 ~ 140 fs. Since 
the hyperfine splittings of |5S 1 / 2 ) and |5 Pi/ 2 ) states (up 
to 3GHz) are much smaller than the laser bandwidth of 
THz order, the hyperfine states can be reduced to a two- 
state system through Morris-shore transform | 2 Tj , giving 



FIG. 3: (Color online) Schematics of the experimental setup. 
Femto-second laser pulses were shaped by an acousto-optic 
programmable dispersive filter (AOPDF) and a 4/-geometry 
zero-dispersion stretcher; and as-produced chirped zero-area 
pulses interacted with cold rubidium atoms ( S5 Rb) in a 
magneto-optical trap. After 3 ns, frequency-doubled laser 
pulses ionized the excited atoms. 


H = 1.46 x 10 " 29 Cm [223. 

Each femto-second laser pulse was divided into two: 
the first (control) pulse after being shaped to a chirped 
zero-area pulse induced the atomic transition, and the 
second (probe) pulse was frequency-doubled to ionize 
atoms in the excited state. The chirped zero-area 
pulse was prepared in two stages: the laser pulse was 
first frequency-chirped, up to c 2 = 6 x 10 4 fs 2 , by 
an acousto-optic programmable dispersive filter (Daz- 
zler from Fastlite) [23] , and the spectrum near the res¬ 
onance was removed by a spectral block in the Fourier 
plane defined by a 4/-geometry Martinez zero-dispersion 
stretcher (homemade) [241125], which consisted of a pair 
of R=500 mm cylindrical mirrors and a pair of gratings 
with 1800 grooves/mm. For the spectral block, we used 
various metal wires of which the spectral width ranges 
from 200 to 1000 GHz in FWHM (or Aw 2 = 7.7 x 10 11 ~ 
3.9 x 10 12 rad/s in Gaussian width). After the two-stage 
of pulse shaping, the laser pulse energy was up to 20 /rJ. 

We used a conventional magneto-optical trap (MOT) 
to spatially isolate the rubidium atoms [26]. By adjust¬ 
ing the diameters of the cooling and repumping laser 
beams, an atomic cloud of about 200 /jm diameter and 
6 x 10 9 cm 3 atom density was prepared. The atoms were 
tightly confined in particular to achieve a uniform laser- 
atom interaction m ■ When the laser pulses were fo¬ 
cused on to the atomic vapor, the laser beam diameter 
of 600 ^m was about 3 times bigger than the diameter of 
the atom cloud. With this diameter ratio, we achieved 
95% of high fidelity for a 7 r-area transform-limited pulse 
excitation. The laser pulse energy (up to 20 /rJ) with 
the given beam diameter was equivalent to 0 ma x = 3 . 57 T. 
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After being interacted with the control pulse (the chirped 
zero-area pulse), atoms in the excited state were ionized 
by the probe pulse, and resulting ions were measured by 
a multichannel plate detector. The overall experimen¬ 
tal cycle controlled by mechanical shutters in laser beam 
lines was 2 Hz to grant the restoration of the MOT. 


V. EXPERIMENTAL RESULTS AND 
DISCUSSION 


Figure [4]( a) shows the experimental result. By varying 
the power and the chirp parameter C 2 of the chirped zero- 
area pulse, the excited-state population was measured. 
The spectral block of A/(FWHM) = 500 GHz (Au >2 = 
1.9x 10 12 rad/s) around the resonance was removed in the 
4/-geometry stretcher. In comparison, the correspond¬ 
ing numerical calculation is shown in Fig. §b). The 
numerical calculation is based on Eq. ©. Since the 
experiment was performed under the condition of non- 
uniform spatial profiles of the laser pulse and the atom 
cloud, the calculation took into account the spatial aver¬ 
aging effect [27], where B^bi i n the x-axis denotes the 
maximum of the Gaussian distribution of 0R a bi- The 
measured populations were calibrated to probabilities by 
using the first peak of Rabi oscillations as a reference. 
Within the available measurement region, the experimen¬ 
tal result shows that the dark CPR region is surrounded 
by the bright CPI region, in a good agreement with the 
theoretical calculation. In addition, we probed the ef¬ 
fect of the width of the spectral hole. The frequency 
chirp was fixed at C 2 = 5 x 10 4 fs 2 , and the spectral hole 
was varied from 200 to 1000 GHz (FWHM). In Fig. |4|c), 
the measured pulse energies (shown with squares) for the 
first CPI peaks are displayed in comparison with the cor¬ 
responding numerical calculation (dots). 

We now consider possible applications of the chirped 
zero-area excitation implicated from the results obtained 
in this study. Since the two-state system interacting with 
the chirped zero-area pulse can undergo, amid the adia¬ 
batic evolution (z-rotation) by the main pulse, Rabi-like 
oscillation due to the non-adiabatic interaction by the 
hole pulse, an interplay between Rabi-like oscillations 
and adiabatic evolutions can be made by shaping the 
spectral width and position of the hole. This approach 
can be a powerful and alternative control means in selec¬ 
tive excitation of, in particular, multi-state systems. In 
a U-type system [2S] of 5S!/ 2 , 5 Pi/ 2 , and 5P 3 / 2 in ru¬ 
bidium, for example, a numerical simulation (not shown) 
with chirped zero-area pulses results in over 97% popu¬ 
lation of the system driven to either excited states, by 
simply changing the pulse intensity only. Also, this con¬ 
trol method implicates that the laser spatial profile is 
also useful for position-dependent selective excitations, 
which can be applied to, for example, atom or ion qubits 
in spatial arrangements mm- 


(a) 
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t/3 
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Spectral block width (THz) 

FIG. 4: (Color online) (a) Experimental result and (b) nu¬ 
merical simulation of the chirped zero-area pulse excitation 
of a cold rubidium atom ensemble. The transition probabili¬ 
ties were plotted as a function of C 2 (linear chirp) and 
Interaction inhomogeneity due to the Gaussian spatial profile 
of the atom ensemble was taken into account for the calcu¬ 
lation of (b). (c) Experimental result of spectral hole-width 
scan. Measured pulse energies for the first CPI (e.g., the stars 
marked in (a) and (b)) are plotted for various spectral hole 
widths (red squares), in comparison with the calculation re¬ 
sult (black dots). 


VI. CONCLUSION 


We have presented our theoretical and experimental 
investigation of the two-state system dynamics under 
the interaction of chirped zero-area pulses. In experi¬ 
ments, we used femto-second laser pulses of which the 























resonant spectral component was removed and then the 
pulses were spectrally chirped. These chirped zero-area 
pulses drove coherent excitations of the resonant 5S]/2- 
5 Pi/ 2 transition of rubidium atoms that were tightly con¬ 
fined in a MOT within a fraction of the laser spatial 
profile to ensure uniform laser-atom interactions. The 
interplay between the adiabatic evolution and Rabi-like 
oscillations, both of which were simultaneously induced 
by the chirped zero-area pulses, has been probed. We 
have shown that the given dynamics can be modeled 
to the Ramsey-type three-pulse interaction, in a good 
agreement with the experimental results. The underly¬ 
ing mechanism behind the observed coherent dynamics 
can be understood based on the nature of the chirped 
zero-area pulse: The chirped zero-area pulse is a sum 
of two pulses with different bandwidtlis, and the broad¬ 
band pulse (the main pulse) makes adiabatic evolution 
even by a small amount of chirp. However, the narrow¬ 


band pulse (the hole pulse) makes either adiabatic or 
non-adiabatic evolution (thus Rabi oscillation). As a re¬ 
sult, the latter switches on or off a Rabi transition amid 
the adiabatic evolution induced by the former. The re¬ 
sult suggests a new design scheme of laser pulse shaping 
towards selective excitation in multi-state atoms, by em¬ 
bedding a Rabi-like oscillation in adiabatic evolutions. 
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